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摘要 :【 目 的] 中华 蜜蜂 Apis cerana cerana 是 一 种 兼 有 生态 价值 和 经 济 价值 的 授粉 昆虫 。 本 研究 拟 
揭示 秦 巴 山区 中 华 蜜蜂 种 群 遗 传 多 态 性 现状 ,探讨 中 华 蜜蜂 的 种 群 分 化 机 制 及 其 影响 因素 。[【 方 
法 】 使 用 8 个 微 卫星 DNA 标记 评估 秦 巴 山区 17 个 样 点 共 979 个 蜂 群 的 中 华 蜜 蜂 遗 传 多 态 性 ,以 长 
白山 中 华 蜜蜂 和 阿坝 中 华 蜜蜂 作为 外 群 进 行 种 群 让 传 分 化 分 析 。[【 结果 】] 秦 巴山 区 中 华 蜜蜂 9596 
的 差异 来 源 于 样 点 内 , 样 点 间 遗 传 分 化 系数 ( FS) 为 0.002 ~0.037, 基 因 流 参数 Nm 为 6.51 ~124.75， 
与 外 群 的 遗传 分 化 分 析 中 均 显 示 秦 巴山 区 中 华 蜜蜂 不 存在 遗传 分 化 。 平 均 期 望 杂 合 度 ( He) 为 
0.6877 x0. 1098 ,平均 观察 杂 合 度 (Ho) 为 0.6364 +0.1367, 平 均 有 效 等 位 基因 数 (Ne) 为 3.7488 圭 
1.6201 个 ,平均 等 位 基因 数 (Na) 为 7.71 +2.52 个 ,多 态 信息 含量 (PIC) 为 0.6418 上 0. 1152 , 香农 
HAI) A 1.5026 + 上 0.3754 ,这 些 参 数 均 表明 微 卫 星 遗 传 多 态 性 丰富 。 中 蜂 圳 状 幼虫 病 流行 过 的 
地 区 ,其 杂 合 度 、 多 态 信 息 含量 以 及 等 位 基因 数 均 显 著 低 于 未 发 病 地 区 。[【 结论 ] 秦 巴山 区 中 华 蜜 
蜂 种 群 数量 大 ,分 布 均匀 ,基因 流水 平 高 ,在 650 km 距离 范围 内 没有 发 生 种 群 分 化 ; 秦 巴 山区 中 华 
蜜蜂 微 卫 星 遗 传 多 态 性 丰富 , 仅 部 分 地 区 中 华 蜜蜂 遗传 多 态 性 受 中 蜂 占 状 幼虫 病 的 选择 压 而 降低 。 
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Population genetic analysis of Apis cerana cerana from the Qinling-Daba 


Mountain Areas based on microsatellite DNA 
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Fuzhou 350002, China) 

Abstract: [ Aim] The Chinese honeybee (Apis cerana cerana) plays vital roles in pollination for agro- 
ecosystem and honey bee products. This study aims to investigate the population differentiation and 
genetic polymorphism among Apis cerana cerana populations across the Qinling-Daba Mountain Areas, 
and to explore possible factors contributing to genetic differentiation and difference of genetic 
polymorphism. [Methods] The genetic relationship was studied, and polymorphism was estimated among 
979 colonies of A. c. cerana from 17 localities combined with two outgroups, A. c. cerana colonies from 
Changbai Mountain and Aba, based on eight microsatellite loci. [ Results] The results showed that up to 
95% of the variation was observed within sampling localities while less variation between localities. The 
genetic differentiation coefficient ( Fst) among 17 localities ranged from 0. 002 to 0. 037, and the gene 
flow (Nm) was from 6. 51 to 124. 75, indicating that no population differentiation was detected. The 
mean expected heterozygosity ( He) , observed heterozygosity ( Ho) , number of alleles ( Na) , number of 


effective alleles ( Ne) , polymorphism information content ( PIC) and Shannon’ s information index (7) 
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were 0.6877 + 0. 1098, 0. 6364 +0. 1367, 7. 71 € 2. 52, 3. 7488 € 1. 6201, 0. 6418 € 0. 1152 and 
1.5026 x0.3754, respectively, suggesting the abundant genetic polymorphism. The values of the genetic 


polymorphism indices including the expected heterozygosity, number of alleles, number of effective 


alleles and polymorphism information content were significantly lower in localities which underwent 


Chinese sacbrood virus pandemic recently than those in the virus uninfected areas. [Conclusion] In the 


context of the large population size and continuous distribution in Qinling-Daba Mountain Areas, there is 


no genetic differentiation among A. c. cerana populations across a distance of 650 km owing to a high 


level of gene flow. And the results also imply that the honeybee in Qinling-Daba Mountain Areas has rich 


genetic polymorphism. 


Significantly decreased polymorphism occurs in some localities that have 


undergone selective pressure from Chinese sacbrood virus prevalence. 


Key words: Apis cerana cerana; genetic 


microsatellite markers 


种 群 遗 传 学 是 研究 种 群 遗 传 分 化 .遗传 结构 LR 
传 特征 .遗传 多 态 性 的 科学 ( 王 崇 云 , 2007) 。 研 究 
种 群 遗 传 分 化 与 特征 , 了解 遗传 资源 现状 是 资源 开 
发 利用 和 遗传 育种 的 前 提 。 遗 传 多 态 性 丰富 是 进化 
的 基础 ,有 利于 生物 适应 新 环境 的 改变 。 评 估 遗 传 
多 态 性 水 平和 探究 其 影响 因素 ,是 保护 和 利用 遗传 
资源 的 必要 条 件 ( 陈 灵芝 和 马克 平 , 2001), 

微 卫 星 标记 广泛 存在 于 真 核 生 物 基因 组 中 ,具有 
共 显 性 多 态 性 高 等 特点 ,被 广泛 地 应 用 于 种 群 遗 传 多 
态 性 的 检测 ( 韩 海底 等 , 2013; Loucifayad et al., 2015) 、 
评估 基因 渗透 ( Harpur et al., 2012, 2013; De la Rúa et 
al., 2013) .人 研究 种 群 遗 传 结构 与 特征 ( Llewellyn et al., 
2008; Oleksa et al., 2011) 分析 遗传 分 化 及 其 机 制 等 
( 张 志 春 等 , 2013; IKIE KSF, 2015) 。 

蜜蜂 是 重要 的 授粉 昆虫 ,对 生物 多 样 性 及 促进 
农 牧 业 增 产 具 有 重要 的 价值 和 意义 (Kremen et al.， 
2007)。 我 国 中 华 蜜蜂 长 期 适应 于 本 土生 态 环 境 ， 
相 比 于 引进 的 西方 塞 蜂 具 有 抗 螨 害 、 耐 寒 善于 利用 
FERTA E KAKER, 2000)。 借 鉴 
西方 蜜蜂 分 子 标记 人 研究 我 国 中 华 蜜蜂 种 群 遗 传 分 
化 ,明确 了 海南 ( 徐 新 建 等 , 2013a) 和 长 白山 ( 张 春 
EE, 2011; TIRE SE, 2013) 的 中 华 蜜蜂 遗传 分 
化 ,探讨 了 环境 异 质 性 以 及 地 理 阻 碍 对 基因 流 的 影 
响 。 然 而 ,地 理 距 离 对 连续 分 布 的 中 华 蜜蜂 种 群 遗 
传 分 化 的 影响 研究 还 未 涉及 。 中 华 蜜蜂 微 卫星 多 态 
性 的 研究 尚 处 于 初始 阶段 , 仅 分 析 了 黄土 高 原 ( 徐 
新 建 等 , 2013b) .沂蒙 山区 ( 陈 唱 等 ,2008a) 南昌 
( 陈 晶 等 , 2008b) fpi pj tL DC (BISE, 2008) RE 
山 ( 朱 翔 杰 等 , 2011a) WERAK, 2011b) \ 海 
南 ( 徐 新 建 等 , 2013a) WE 曹 联 飞 等 , 2013) 和 长 
白山 (于 涉 龙 等 , 2013 ) 等 局 部 地 区 ,中 华 蜜蜂 微 卫 



































出 









































differentiation ; 


genetic polymorphism; population; 


星 的 期 望 杂 合 度 水 平 为 0.21 ~0.73。 外 来 引种 、 疾 
病 ( 于 泳 龙 等 , 2013; 徐 新 建 等 , 2013a) 、 蜜 粉 源 ( 陈 
唱 等 , 2008a, 2008b; 徐 新 建 等 , 2013b) 、 人 工 选 育 
( 朱 翔 杰 等 , 2011b; 曹 联 飞 等 , 2013 ) 等 因素 可 能 影 
响 中 华 蜜 蜂 的 遗传 多 态 性 。 

秦 巴 山区 位 于 我 国 的 中 西部 , 介 于 102°54' ~ 
112?40' E 和 30° 50' ~ 34°59' N 之 间 , 总 面积 约 
22.23 x 10* km ,地 跨 陕 西 .甘肃 四川、 湖北 、 河 南 及 
重庆 等 省 市 ( 白 红 英 , 2014)。 秦 巴山 区 植被 良好 ,中 
华 蜜蜂 数量 较 多 ,是 研究 大 面积 相似 生境 条 件 下 ,中 
华 蜜蜂 种 群 遗 传 多 态 性 和 遗传 分 化 的 理想 场所 。 

本 研究 采用 微 卫 星 标记 ,分 析 秦 巴山 区 中 华 蜜 
蜂 的 种 群 分 化 与 遗传 结构 ,尝试 解决 在 生境 相似 、 中 
华 蜜蜂 连续 分 布 的 条 件 下 ,地 理 距离 与 种 群 分 化 的 
关系 。 评 估 秦 巴山 区 中 华 蜜蜂 微 卫 星 多 态 性 水 平 ,分 
析 中 蜂 嘻 状 幼虫 病 对 中 华 蜜蜂 遗传 多 态 性 的 影响 ,以 
便 全 面 擎 握 秦 巴山 区 中 华 蜜蜂 遗传 资源 现状 与 特征 ， 
为 利用 和 保护 中 华 蜜 蜂 遗 传 资源 奠定 基础 理论 。 









































1 材料 方法 


1.1 供 试 昆 虫 

在 秦 巴 山区 17 个 采样 地 点 , 共 979 个 原始 饲养 
的 中 华 蜜蜂 蜂 群 中 采集 工蜂 样本 。 原 始 饲 养 的 中 华 
蜜蜂 处 于 半 野 生 状 态 , 不 进行 人 工 育 王 , 相 比 活 框 饲 
养 的 中 华 蜜蜂 样本 的 代表 性 好 。 采 样 时 ,打开 原始 
蜂 集 在 梨 内 的 边缘 取样 或 在 集 门口 提取 工蜂 。 秦 岭 
山脉 中 的 陕西 西南 部 和 甘肃 东南 部 的 中 华 蜜 蜂 在 采 
样 前 5 年 内 曾 发 生 严重 的 中 蜂 吉 状 幼虫 病 ( 表 1) 。 
以 长 白山 和 阿坝 各 1 个 采样 地 点 的 样本 作为 研究 遗 
传 分 化 的 外 群 对 照 。 




















































































































3 期 郭 慧 萍 等 : 秦 巴 山区 中 华 蜜蜂 种 群 微 卫星 DNA 遗传 分 析 339 
表 1 本 研究 中 华 蜜蜂 19 个 采样 地 点 的 样本 信息 
Table 1 Information of total nineteen collecting sites of Apis cerana cerana in this study 

样本 代码 采集 地 点 采样 群 数 经 纬度 Id (m) 

Sample code Collecting site Number of colonies sampled Latitude and longitude Altitude 
MJ 甘肃 麦 积 Maiji, Gansu * 63 34*17.601'N; 106?01.447'E 1 610 
HX 甘肃 微 县 Huixian, Gansu * 50 33559. 793'N; 106^05.034'E 1 167 
FX 陕西 风 县 Fengxian, Shaanxi * 52 34*07.977'N; 106258. 871'E 1 641 
WB 陕西 渭 滨 Weibin, Shaanxi * 57 34°20. 536'N; 107°01. 447'E 752 
TB 陕西 太白 Taibai, Shaanxi " 54 34°03. 603'N; 107°23. 983'E 1 655 
LB 陕西 留 坝 Liuba，Shaanxi” 53 33?34.557'N; 106°57. 079'E 957 
NS 陕西 宁 陕 Ningshan, Shaanxi * 55 33°40. 580'N; 108?18.838'E 1 390 
NQ 陕西 宁 强 Ningqiang, Shaanxi * 59 32*47.229'N; 10610. 816'E 868 
ZP 陕西 镇 平 Zhenping, Shaanxi 56 32?05.404'N; 109°25. 006'E 1 130 
SX 河南 陕 县 Shanxian, Henan 60 34?31.354'N; 111°22. 460'E 782 
NY 河南 南阳 Nanyang, Henan 77 33?16.311'N; 111?31.056'E 224 
WY 四 川 万 源 Wanyuan, Sichuan 65 31°49. 880'N; 108?21.326'E 1 529 
QC 四 川 青川 Qingchuan, Sichuan 63 32°35. 554'N; 104?44. 134'E 1 535 
NJ 四 川南 江 Nanjiang, Sichuan 62 32?37.226'N; 107°00. 584'E 976 
WS 庆 巫 山 Wushan, Chongqing 52 31?12.620'N; 109?45.341'E 691 
CK 庆 城 口 Chengkou, Chongqing 50 31?40. 073'N; 108°11. 000'E 881 
YY 庆 云 阳 Yunyang, Chongqing 51 30°43.705'N; 108?47. 534'E 543 
CBS 吉林 长 白山 Changbai Mountain, Jilin 38 42*32.581'N; 128930. 563'E 725 
AB 四 川 阿 坝 Aba, Sichuan 67 31?53. 271'N; 102°07. 934'E 2717 





* AK Z5 Jp EP ESI E Having suffered from Chinese sacbrood virus. 


1.2 基因 组 DNA 的 提取 与 扩 增 


每 





使 用 1 头 工蜂 样本 进行 基因 组 的 提取 与 微 


卫星 扩 增 (Harpur et al., 2012; Chahbar et al., 
2013) 。 利 用 Ezup 柱 式 基因 组 DNA 抽 提 试剂 盒 | 生 
工 生物 工程 (上 海 ) 股份 有 限 公 司 ] 提取 蜜蜂 基因 
组 ,利用 ExTaq Hot Start Version 扩 增 试剂 盒 [ 宝生 








物 工程 (大 连 








研究 , PCR 产物 经 ABI 3730xl 自动 测序 仪 (美国 








) 有 限 公 司 , 大 连 ] 对 8 个 微 卫 星 标记 
3 3$ ( Michel et al., 2007; Takahashi et al., 2009) 
(3€ 2) , PCR 的 扩 增 体系 参照 徐 新 建 等 (2013a) 的 





上 


ABI) 进行 毛细 管 电泳 和 等 位 基因 分 型 ( CeneMapper 


4.0)。 


A2 8 对 微 卫 星 标记 引物 信息 


Table 2 Primers of eight microsatellite markers 





位 点 引物 序列 (5 -37) 退火 温度 (%C ) GenBank 登录 号 
Locus Primer sequences Annealing temperature GenBank accession no. 
an ee Tcz— uo 
ASH F: GCTACACTACCGCGACCTGCA 58 AB363803 
R: TACGCTCCCTTACTCCCGCTG 
he2 ,GAAGACAGAGTCAAAACCGCA s AB363804 
Àc5 F: TGGCCGAAGAACGGTGTTACG 58 AB363806 
R: GAAAGAAACGAGTCGCTGGTGT 
F: AACCTICTTCGCCACCTCCAA 
AER R: GTCTGAACGAAAGAAACAGCA 38 AB363809 
Ac35 F: TTCCAGCAGGAAGTGACGGTG 58 AB363814 
R: CGGGAAACTCGTCATTTTCGA 
Ap085 F: GATCAAACACACAAACGAAAGC 58 AJ509359 
R: ACCGGAAGCCTAATCAAGG 
ATioi 4 GCGTTCCAAGTGAATGAACA 58 AJ509549 


: GTTGGCTATTTTCGTATCGC 
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1.3 数据 统计 与 分 析 

利用 Excel Micro-Satellite Toolkit 3. 1 ( Park, 
2001) 软件 计算 观察 杂 合 度 ( Ho) BEER BE ( He) 
和 等 位 基因 数 (Na) ,PopGene 1.31 软件 分 析 样 点 有 
效 等 位 基因 数 ( we) 和 香农 指数 (1)。GenAlEx 6.5 
软件 进行 主 成 分 分 析 (PCA ) 和 AMOVA 分 析 
(Peakall and Smouse, 2012) , Genepop on Web 在 线 
软件 分 析 样 点 间 的 遗传 分 化 系数 (Fsi) ( Rousset, 
2008) ,由 遗传 分 化 系数 Fst 计算 基因 流 参数 Nm, 
Nm=1/4(1/Fst -1) (Slatkin et al., 1989; Wang et 
al, 2015), fi Hj Poptree 2 构建 系统 发 生 树 ， 
Structure 2. 3. 4 软件 分 析 基 因 遗 传 结构 (Pritchard et 
al., 2000) 。 


2 结果 


2.1 秦 巴 山区 中 华 蜜蜂 种 群 遗传 分 化 
根据 秦 巴 山区 中 华 蜜蜂 样 点 间 的 遗传 分 化 系 
数 、 遗 传 距离 遗传 结构 分 析 、 主 成 分 分 析 和 系统 发 
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生 树 构建 等 研究 结果 ,证 明 各 样 点 间 的 中 华 蜜蜂 均 
未 发 生 种 群 遗 传 分 化 。 

秦 巴 山区 中 华 蜜蜂 样 点 间 遗 传 分 化 系数 为 
0. 002 ~ 0. 037 ; 遗传 距离 为 0.029 ~0.215( 表 3)。 
AMOVA 分 析 结 果 表 明 高 达 95% 的 变异 来 源 来 自 于 
样 点 内 部 个 体 间 。 根 据 遗 传 分 化 系数 计算 的 基因 流 
参数 Nm 值 为 6.51 ~ 124. 75 ,表明 各 样 点 间 的 基因 
流 大 ,没有 导致 遗传 分 化 。 

将 吉林 省 长 白山 和 四 川 省 阿坝 样本 作为 外 群 与 
秦 巴 山区 中 华 蜜蜂 样 点 进行 遗传 结构 分 析 , 秦 巴山 
区 各 样 点 的 中 华 蜜蜂 遗传 结构 相似 ,与 2 个 外 群 的 
遗传 结构 差别 明显 (图 1)。 

主 成 分 分 析 表 明 秦 巴山 区 的 中 华 蜜蜂 没有 分 
化 。 第 1 主 成 分 与 第 2 主 成 分 贡献 率 分 别 为 
57.44% 和 27.96% , 秦 巴 山区 中 华 蜜蜂 紧密 的 聚 在 
一 起 ,与 长 白山 、 阿 坝 中 华 蜜蜂 分 化 明显 (图 2)。 

聚 类 分 析 表 明 秦 巴山 区 的 中 华 蜜蜂 没有 分 化 ， 
在 聚 类 树 上 聚 在 一 起 ,与 长 白山 阿坝 中 华 蜜蜂 明显 
分 化 (图 3)。 
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图 1 秦 巴 山区 \ 长 白山 和 阿坝 中 华 蜜蜂 的 遗传 结构 分 析 


Fig. 1 





Genetic structure analysis of Apis cerana cerana from the Qinling-Daba Mountain Areas, Changbai Mountain and Aba 


样本 代码 同 表 1。The sample codes same as in Table 1. 





0 阿坝 种 群 Apa population 











A 长 白山 种 群 Changbai Mountain 
population 

O 秦 巴 山区 种 群 Populations from A 

-0.04 Qinling - Daba Mountain Areas 








0.08 
-0.15  -0.10 -0.05 0 0.05 0.10 0.15 
主 成 分 1 
Factor 1 
图 2 秦 巴 山区 ,长 白山 和 阿坝 中 华 蜜蜂 主 成 分 分 析 


Fig. 2 Principal component analysis of Apis cerana cerana 








from the Qinling-Daba Mountain Areas, Changbai 


Mountain and Aba 


2.2 秦 巴 山区 中 华 蜜蜂 的 多 态 性 分 析 

秦 巴 山区 中 华 蜜蜂 从 期 望 杂 合 度 .观察 杂 合 度 、 
等 位 基因 数 、 有 效 等 位 基因 数 、 多 态 信息 含量 和 香农 
指数 等 6 个 多 态 性 参数 分 析 , 微 卫星 遗传 多 态 性 水 
平 较 高 。 严 重病 害 能 够 降低 中 华 蜜蜂 的 遗传 多 
态 性 。 
2.2.1 秦 巴 山区 中 华 蜜蜂 的 多 态 性 : 秦 巴 山区 中 华 
蜜蜂 期 望 杂 合 度 为 0.6497 ~0.7288 ,观察 杂 合 度 为 
0.5686 ~0.7078 ,等 位 基因 数 为 6.75 ~ 8.75 ,有效 
等 位 基因 数 为 3. 2659 ~ 4. 4271, 多 态 信息 含量 为 
0. 6012 ~ 0. 6791, 香农 指数 为 1. 3766 ~ 1. 6813 
( 表 4)。 这 6 项 多 态 性 参数 均 能 表明 秦 巴 山区 中 华 
蜜蜂 的 多 态 性 丰富 。 
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图 3 秦 巴 山区 \ 长 白山 和 阿坝 中 华 蜜蜂 UPGMA RKI 
Fig. 3 UPGMA dendrogram of Apis cerana cerana in the Qinling-Daba Mountain Areas, Changbai Mountain and Aba 
标尺 表示 遗传 距离 。The scale bar represents the genetic distance. 
表 3 秦 巴 山区 样 点 间 中 华 蜜蜂 种 群 遗传 距离 和 遗传 分 化 系数 
Table 3 Genetic distance and genetic differentiation coefficient between samples of Apis cerana cerana 
in Qinling-Daba Mountain Areas 
0 NS ZP FX TB WB LB NQ MJ HX SX NY NJ WY QC CK WS YY 
NS - 0.042 0.046 0.058 0.038 0.163 0.041 0.048 0.043 0.067 0.042 0.043 0.069 0.084 0.045 0.052 0.058 
ZP 0.013 - 0.062 0.062 0.041 0.165 0.047 0.063 0.054 0.052 0.044 0.035 0.091 0.080 0.045 0.059 0.049 
FX 0.014 0.023 - 0.061 0.034 0.166 0.052 0.060 0.062 0.081 0.071 0.057 0.084 0.091 0.056 0.071 0.073 
TB 0.023 0.023 0.022 - 0.057 0.166 0.059 0.069 0.072 0.080 0.063 0.060 0.088 0.100 0.075 0.079 0.088 
WB 0.013 0.013 0.006 0.028 - 0.161 0.040 0.051 0.049 0.067 0.053 0.046 0.080 0.090 0.050 0.061 0.051 
LB 0.004 0.013 0.013 0.010 0.013 - 0.161 0.167 0.166 0.182 0.157 0.163 0.215 0.169 0.172 0.173 0.174 
NQ 0.010 0.018 0.016 0.021 0.011 0.008 - 0.052 0.048 0.066 0.040 0.050 0.082 0.075 0.050 0.057 0.065 
MJ 0.012 0.026 0.019 0.019 0.016 0.005 0.010 - 0.029 0.075 0.051 0.060 0.096 0.081 0.059 0.068 0.069 
HX 0.006 0.015 0.021 0.020 0.018 0.004 0.007 0.005 - 0.063 0.047 0.059 0.105 0.075 0.057 0.061 0.066 
SX 0.019 0.012 0.021 0.027 0.023 0.012 0.025 0.028 0.023 - 0.057 0.060 0.130 0.098 0.067 0.082 0.087 
NY 0.017 0.010 0.026 0.019 0.023 0.008 0.011 0.015 0.014 0.016 - 0.041 0.093 0.067 0.045 0.044 0.048 
NJ 0.017 0.003 0.021 0.019 0.023 0.015 0.028 0.030 0.022 0.019 0.014 - 0.087 0.078 0.050 0.040 0.046 
WY 0.027 0.019 0.032 0.029 0.032 0.030 0.024 0.037 0.027 0.033 0.021 0.020 — 0.134 0.080 0.093 0.100 
QC 0.023 0.002 0.029 0.029 0.019 0.014 0.020 0.028 0.019 0.021 0.010 0.013 0.026 - 0.067 0.077 0.076 
CK 0.015 0.016 0.018 0.022 0.015 0.006 0.007 0.013 0.005 0.025 0.014 0.020 0.020 0.007 - 0.062 0.054 
WS 0.010 0.005 0.019 0.026 0.015 0.010 0.019 0.022 0.012 0.015 0.009 0.005 0.016 0.008 0.009 - 0.050 
YY 0.015 0.007 0.026 0.030 0.015 0.011 0.023 0.025 0.015 0.030 0.012 0.009 0.016 0.004 0.009 0.002 Z 








对 角 线 上 三 角 为 遗传 距离 (De), 下 三 








below the diagonal. 
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对 多 态 性 的 影响 : 
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为 遗传 分 化 系数 (Fst) o Genetic distance is above the diagonal, while the genetic differentiation coefficient is 


(t = -4. 50899, P =0. 00042) 、 有 效 等 位 基因 数 
(t= -3.07693, 已 =0.00767) 极 显著 地 低 于 未 发 病 
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显著 低 于 未 发 病 样 点 。 
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XA 秦 巴 山区 中 华 蜜蜂 遗传 多 态 性 
Table 4 Genetic TT of Apis cerana cerana from the Qinling-Daba Mountain Areas 
样本 代码 期 望 杂 合 度 见 察 杂 合 度 多 态 信息 合 量 。 等 位 基因 数 。 有效 等 位 基因 数 FRA 

Sample code He Ho PIC Na Ne I 
NS 0.6618 +0.1411 0.6383 +0.1119 0.6202 +0. 1426 7.25 +2.38 3.8483 + 2.0735 1.5062 +0. 3759 
FX 0.6721 x0.1189 0.6168 x0. 1098 0.6216 x0. 1258 6.75 x2.25 3.5246 + 1. 6467 1.4074 +0. 3598 
TB 0.6567 x 0.1087 0.6196 x 0.0785 0.6095 x0. 1097 7.50 x2.14 3.4561 + 1.3654 1.4787 x0.2726 
WB 0.7010 x0. 1066 0.6549 x0. 1341 0.6567 x0.1132 7.50 +2.88 3.9319 «1.3680 1.5768 +0. 3595 
LB 0.6759 +0. 1144 0.7078 +0. 1300 0.6363 +0. 1194 8.29 +3.04 3.4621 +2. 0440 1.3766 +0. 6641 
NQ 0.6620 +0. 1051 0.6114 +0. 1298 0.6203 +0. 1097 7.75 +2.49 3.5964 + 1.5375 1.5385 +0. 3132 
MJ 0.6547 +0. 1242 0.5686 +0. 1626 0.6107 +0. 1252 7.25 +1.49 3.2756 +1.2817 1.3954 +0. 3157 
HX 0.6497 +0. 1214 0.5861 +0. 1926 0.6012 +0. 1259 7.13 x1. 73 3.2659 + 1.4573 1.3774 x0. 3227 
ZP 0.7143 +0. 0847 0.6959 +0. 1444 0.6691 +0. 0916 8.50 +2.67 3.9763 + 1.4265 1.6060 +0. 2992 
SX 0.7207 x 0.0791 0.6212 x0. 1448 0.6791 x0. 0863 8.75 +2.38 4.1691 x 1.1746 1.6813 +0. 2544 
NY 0.7028 +0. 0934 0.6810 +0. 1281 0.6587 +0. 1017 7.88 +2.80 3.8312 + 1.7262 1.5204 +0. 3272 
NJ 0.6972 +0. 1165 0. 6023 +0. 1081 0.6538 +0. 1188 8.00 +2.51 3.7509 + 1.4690 1.5119 +0. 3438 
WY 0.6755 x0. 1383 0.5948 x0. 2052 0.6263 +0. 1528 7.75 +3.45 3.8500 +2. 0319 1.4386 +0. 3438 
QC 0.7288 x 0.0932 0.6465 x 0.0975 0.6572 x0. 1001 8.00 +3.55 3.5500 x 1. 8817 1.3778 x0. 5896 
CK 0.6862 x0.1110 0.6704 x0. 1253 0.6446 x 0.1186 7.50 x 2.93 3.77034 € 1. 7841 1.5021 +0. 3585 
WS 0.7196 x 0.0962 0.6915 +0. 1954 0.6766 +0. 0987 7.38 +1.41 4.1107 + 1.5431 1.5881 +0. 3073 
YY 0.7127 x0.1146 0.6109 +0. 1258 0.6693 x0.1188 7.88 x2. 70 4.4271 x1. 7312 1.6611 x0. 3383 

平均 值 Mean 0.6877 +0. 1098 0.6364 +0. 1367 0.6418 +0. 1152 Tatl 322:52 3.7488 + 1. 6201 1.5026 +0. 3754 

表 中 数据 为 平均 值 + 标准 差 。Data in the table are means + SD. 

发 生病 害 的 样 点 , 微 卫 星 8 个 位 点 等 位 基因 均 ”比较 单个 位 点 的 6 个 多 态 性 参数 ,在 Ac-2 位 点 均 显 





发 生 不 同 程度 丢失 ,其 中 Ac2 和 ATIOL 等 位 基因 
丢失 最 多 ,分 别 丢失 50% 和 35.2% 。 发 病 组 与 未 发 
病 组 多 态 性 6 个 多 态 性 参数 均值 比较 均 显著 降低 ; 


表 5 


中 蜂 圳 状 幼虫 病 发 病 区 与 未 发 病 区 中 华 蜜蜂 





著 低 于 未 发 病 区 ， Ac-2 和 Ac-5 位 点 的 多 态 RAS 
量 及 Ac-2, Ap085, AT101 和 Ac-1 位 点 观察 杂 合 度 


车 低 于 未 发 病 区 ( 表 5) 。 
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Table 5 Differences in genetic polymorphism of eight microsatellite loci in Apis cerana cerana between 


Chinese sacbrood disease infected and uninfected areas 
























































等 位 基因 种 类 期 望 杂 合 度 观察 杂 合 度 多 态 信息 含 等 位 基因 数 有效 等 位 基因 数 香农 指数 
位 点 All" He Ho PIC Na Ne I 
Locus — A “未 发 病 ”发 病 “未 发 病 ”发病 “未 发 病 ”发病 “未 发 病 ”发 病 ” 未 发 病 ”发病 “未 发 病 ”发病 ”未 发 病 
Infected Uninfected Infected Uninfected Infected Uninfected Infected Uninfected Infected Uninfected Infected Uninfected Infected Uninfected 
Ac-35 12 13 0.77 0.79 0.77 0.73 0.72 0.75 9.00 9.89 4.63 4.85 1.72 1.80 
Ac-2 9 18 0.57* 0.71” 0.57* 0.70” 0.53" 0.67™ 7.00 8.33 2.46" 3.61" 1.25"  1.55* 
Ap085 15 18 0.85 0.86 0.74* 0.80” 0.82 0.83 11.50 13.00 6.77 7.23 2.11 2.17 
Ac-26 10 13 0.69 0.69 0.65 0.69 0.64 0.64 7.13 6.89 3.30 3.52 1.45 1.39 
Ac-5 10 11 0.55 0.62 0.57 0. 60 0.50* 0.58* 6.88 7.33 2.49 2.94 1.20 1.36 
A313 17 18 0.58 0.59 0.58 0.56 0.55 0.55 10.38 10.11 2.58 2.60 1.42 1.38 
ATIOI 11 17 0.72 0.75 0.42* 0.44* 0.67 0.67 7.38 8.44 3.37 3.70 1.33 1.45 
Ac-1 6 7 0.62 0.65 0.65* 0.64* 0.55 0.59 9:25 5.44 2.76 2.99 1.19 1.25 
Ton 11.25 14.38 0.67" 0.71" 0.62* 0.65* 0.62" 0.66" $8.06" 58.68" 3.55" 3.93" 1.46 1.54 
"P«0.05; "P«0.0l. 
大 ,会 由 于 基因 流 的 逐渐 减少 而 导致 遗传 分 化 
+ * 一 | ——— 
3 讨论 (Barton et al., 2010) ,对 于 不 同 生物 而 言 多 远 的 距 


3.1 地 理 距离 与 中 华 蜜蜂 种 群 遗传 分 化 的 关系 
理论 上 ,大 范围 连续 分 布 的 种 群 ,随地 理 距 离 加 





J 


离 会 产生 遗传 分 化 仍 处 于 探索 阶段 。 随 分 子 数据 分 
析 地 理 隔 离 方 法 的 出 现 (Bohanak ，2002 ) ,该 领域 的 
研究 受到 较 多 关注 ,特别 是 营 固 着 生活 
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的 生物 (Hurtrez-Bousses et al., 2010; Husemann et 
al., 2012; Gu et al., 2015 )。 秦 巴山 区 蜜 粉 源 植物 
丰富 ,主要 以 中 山 和 低 山 为 主 ( 白 红 英 , 2014) ,生境 
相似 ,地 理 距 离 是 唯一 有 可 能 影响 中 华 蜜蜂 种 群 遗 
传 分 化 的 因素 ,因此 ,是 研究 遗传 分 化 与 地 理 距 离 的 
天 然 理 想 实 验 场所 。 

Wright(1978) 提出 , 当 Fst <0. 05 ,表明 种 群 间 
弱 分 化 或 没 分 化 。 本 研究 中 ,所 有 样 点 间 遗 传 分 化 
系数 Fst 值 为 0.002 ~ 0.037 ,直线 距离 最 远 的 四 川 
青川 和 河南 陕 县 间 约 为 650 km ,两 者 间 Fst 值 为 
0.021, 四 川 万 源 与 甘肃 麦 积 Fsi 值 最 大 ,为 0.037, 
均 小 于 0.05, 因 此 认为 秦 巴 山区 中 华 蜜蜂 未 发 生 遗 
传 分 化 。 同 时 , AMOVA 分 析 显 示 9596 的 变异 来 源 






































0. 1427 ~0.8880, 均 低 于 秦 巴 山区 在 2 个 位 点 上 的 
期 望 杂 合 度 (0.7337 和 0. 8542) ,说明 秦 巴山 区 中 
华 蜜蜂 微 卫星 遗传 多 态 性 丰富 ,是 我 国 中 华 蜜蜂 遗 
传 资源 保护 和 利用 的 重要 地 区 。 
3.3 ”中 蜂 吉 状 幼虫 病 导致 秦 巴 山区 中 华 蜜蜂 遗传 
多 态 性 降低 的 原因 

影响 蜜蜂 微 卫 星 遗 传 多 态 性 的 因素 很 多 ,主要 
有 选择 压力 、 突 变迁 移 扩散 等 (Zayed,2009 ) ,选择 
压力 对 种 群 遗传 多 态 性 的 影响 又 是 生物 学 家 普遍 关 
注 的 问题 (Harpur et al., 2012; Jara et al., 2015 ) 。 
中 蜂 品 状 幼虫 病 是 严重 危害 中 华 蜜蜂 的 病毒 性 疾 
病 , 流 行 时 蜂 群 的 死亡 率 常 达 50% 以 上 ( 王 双 进 ， 
2008 ) 。 这 种 毁灭 性 疾病 的 流行 对 于 中 华 蜜蜂 的 生 









































于 样本 内 , 肾 类 分 析 和 遗传 结构 分 析 结 果 也 支持 了 
这 个 观点 。 基 因 流 参数 Nm > 1 表示 基因 流 足 以 防 
止 遗 传 分 化 的 发 生 ( Slatkin and Barton, 1989; Wang 
et al., 2015) ,在 西方 蜜蜂 的 研究 中 Nm >2 即 表 示 
不 存在 遗传 分 化 (Nikolova et al., 2015) , 而 本 研究 
中 Nm 为 6.51 ~124.75。 说 明 在 种 群 连续 分 布 , 环 
境 均 质 的 条 件 下 ,中 华 蜜蜂 种 群发 生 遗 传 分 化 的 距 
离 要 大 于 650 km。 

中 华 蜜蜂 婚 飞 常 发 生 在 10 km 以 外 ,平原 可 达 
18 km ÆA CJ] vule, 2002), 秦 巴山 区 中 华 蜜蜂 种 
群 数量 大 , 且 连 续 分 布 ,未 产生 分 化 的 主要 原因 很 可 
能 是 样 点 间 存 在 较 高 的 基因 流 。 我 国 中 华 蜜蜂 分 布 
生态 区 域内 基本 不 存在 比 秦 巴山 区 更 大 的 同 质 环 
境 ,我 们 推测 在 中 国 的 生态 环境 下 ,在 中 华 蜜蜂 连续 
分 布 区 域内 仅 由 地 理 距离 因素 导致 的 中 华 蜜蜂 种 群 
遗传 分 化 的 可 能 性 很 小 。 
3.2 秦 巴 山区 中 华 蜜蜂 的 遗传 多 态 性 水 平 

种 内 的 遗传 多 态 性 是 一 个 物种 对 抗 外 界 干扰 能 
否 成 功 反 应 的 决定 因素 ( 陈 灵芝 和 马克 平 , 2001), 
对 于 培育 新 的 优良 品种 具有 重要 意义 。 期 望 杂 合 度 
是 度量 种 群 遗 传 多 态 性 常用 的 参数 , Takezaki 和 Nei 
(1996) 提 出 微 卫星 计算 的 期 望 杂 合 度 在 0.3 ~0.8 
之 间 可 说 明 种 群 遗 传 多 态 性 较 高 。 厅 巴山 区 中 华 蜜 
蜂 微 卫星 遗传 多 态 性 平均 期 望 杂 合 度 为 0.6877 ,在 
已 经 报道 的 文献 中 中 华 蜜蜂 期 望 杂 合 度 为 0.21 ~ 
0.73( 徐 新 建 等 , 2013a, 2013b; FSE, 2013), 
且 本 研究 中 使 用 的 AT101 和 Ap085 位 点 在 长 白山 
(TAES, 2013) 武夷山 ( 朱 翔 杰 等 , 2011a) \ 海 
南 ( 徐 新 建 等 , 2013a) ` 福 建 ( 朱 翔 杰 等 , 2011b) 、 南 
昌 ( 陈 晶 等 , 2008b) 等 地 区 的 研究 中 也 有 涉及 ,2 个 
位 点 的 期 望 杂 合 度 分 别 为 0. 1780 ~ 0. 7001 和 

































































存 发 展 而 言 是 一 种 较为 强烈 的 选择 压 , 疾 病 对 遗传 
多 态 性 的 影响 根本 原因 是 个 体 的 死亡 造成 等 位 基因 
ER ,最 终 导致 遗传 多 态 性 的 降低 。 本 研究 区 域 种 
群 数量 大 ,总 体 遗 传 多 态 性 丰富 ,但 流行 过 中 蜂 者 状 
幼虫 病 的 样 点 , 微 卫 星 遗传 多 态 性 显著 低 于 未 发 病 
区 域 ,部 分 位 点 出 现 大 量 的 等 位 基因 丢失 的 现象 ,如 
Ac-2 在 发 病 区 与 未 发 病 区 等 位 基因 数 由 18 种 减少 
为 9 种 ,充分 说 明 中 蜂 圳 状 幼虫 病 的 选择 强度 足以 
降低 中 华 蜜蜂 的 遗传 多 态 性 。 
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